NASA TECHNICAL 

MEMORANDUM 


NASA TM X- 62,227 


r- 

CM 

CM 

CD 


I- 

< 

CO 

< 

z 


THE INDUCED MAGNETOSPHERE OF THE MOON 
I. THEORY 


G. Schubert 

University of California, Los Angeles, Calif. 90024 
C. P. Sonett 

Ames Research Center, Moffett Field, Calif. 94035 


K. Schwartz 
20858 Collins Street 
Woodland Hills, Calif. 91364 

H. J. Lee 

University of California, San Diego 
La Jolla, Calif. 92037 


January 1973 


(NASA-TM-X-62227) THE INDUCED 
MAGNETOSPHERE OF THE MOON. Is THEORY 
(NASA) 59 p HC $5.00 CSCI 03B 



Unclas 

G3/30 54909 



The Induced Magnetosphere of the Moon 
I. Theory 


G. Schubert 

Department of Planetary and Space Science 
University of California, Los Angeles 
90024 


C.P. Sonett 

National Aeronautics & Space Administration 
Ames Research Center 
Moffett Field, California 
94035 


K. Schwartz 

20858 Collins Street 
Woodland Hills, California 
91364 


H.J. Lee 

Department of Applied Physics & Information Science 
University of California, San Diego 
La Jol la, Cal ifornia 
92037 



ABSTRACT 


An analytic solution for the magnetic field in the space 
defined by a spherical Moon and its downstream cylindrical 
cavity formed by the solar wind is derived for interplanetary 
magnetic fields both parallel and perpendicular to the cavity 
axis. By superposition, the solution is obtained for arbitrary 
orientations of the interplanetary field. The theory is 
quasistatic and is formulated in terms of a sea 1 a r magnet i c 
potential. Thus the Moon model consists of a core of 
arbitrary size and infinite electrical conductivity surrounded 
by a nonconducting shell; the cavity volume is also assumed to 
be nonconducting. The variation of the magnetic field on the 
lunar surface, both on the sunlit and on the dark side hemispheres, 
and on the cavity boundary is presented for various values of 
core radius. The solution also yields the distribution of 
currents on the lunar sunlit surface and the surface of the 
cavity. Theoretical transfer functions are also presented 
and their variations with position on the lunar surface and 
with core size are discussed. The results are compared with 
those of spherically symmetric approximations to lunar magnetic 
i nduct i on . 
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I NTRODUCTION 

Quantitative measurements of lunar el ectromagneti c in- 
duction by the Apollo 12 surface magnetometer have provided 
a considerable stimulus towards understanding the role played 
by the asymmetric plasma environment of the Moon in this 
phenomenon. Magnetic field fluctuation data, obtained when 
the surface magnetometer is on the front side, i .e. the side 
directly exposed to the solar wind or magnetosheath plasma, 
have been interpreted in terms of a spherically symmetric 
theory which assumes no scattered electromagnetic radiation 
[Sonett et al., 1971a, b; Sonett et al . , 1972a]. The analytic 
formulation of this theory was carried out by Schubert and 
Schwartz [1969, 1972], Schwartz and Schubert [1969] and Blank 
and Sill [1969]. The no scattering assumption is equivalent 
to assuming that magnetic fields associated with induced 
currents in the Moon are confined to the lunar volume. On 
the front side this condition is approximately correct since 
waves cannot travel upstream in the supermagnetosoni c plasma. 
Further, magnetometer experiments on the lunar sunlit surface 
have shown that the fields induced in the lunar interior by 
the action of fluctuations in the interplanetary magnetic 
field are strongly confined by the momentum flux of the solar 
wind [Sonett et al., 1971c]. However, because of the inter- 
ruption of the solar wind flow by the Moon and the consequent 
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lack of plasma dynamic pressure on the dark lunar hemisphere, 
it is likely that magnetic fields induced in the lunar interior 
leak into the cavity. 

Magnetic field fluctuation data, obtained when the surface 
magnetometer is on the far side within the diamagnetic cavity, 
have been interpreted in terms of a spherically symmetric theory 
which assumes vacuum scattering [Dyal and Parkin, 1971; Dyal 
et al., 1972], The theoretical formulation of the spherically 
symmetric vacuum scattering problem is a classic exercise in 
electrodynamics (see e.g. Stratton, 1941). However, the in- 
duced magnetic fields are in fact essentially confined to the 
body of the Moon and the space of the cavity. Approximate con- 
finement to the interior of the cavity occurs since downstream 
travelling scattered waves in the supermagnetoson i c solar wind 
are limited to the interior of the Mach cone extending down- 
stream from the limb. A self-consistent confinement current 
layer therefore exists along the cavity walls and distorts that 
part of the lunar induced field found on the dark hemisphere. 

The confinement current system is conceptually distinct from 
the currents associated with the cavity diamagnetism [Colburn 
et al., 1967]. Both the confinement currents and the diamagnetic 
currents are the result of discontinuities in the tangential 
magnetic field. The model discussed in this paper is restricted 
to confinement currents. 

Lunar magnetic induction at hydromagneti c frequencies is a 
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global phenomenon; it follows that a theory for the induction 

\ 

should account for asymmetric plasma conditions surrounding the 
Moon. In this pap^r, Part I, we present an asymmetric induction 
theory and compare the results with those of the two spherically 
symmetric approximations. In Part II, lunar surface magnetometer 
data will be discussed and interpreted using the asymmetric 
theory. We note in this connection that Reisz et al . [1972] have 
considered the asymmetric induction problem for a cylindrical 
Moon. Their solution was obtained using a numerical finite 
difference scheme but the lack of detail in their paper makes 
an evaluation difficult [Sonett et al., 1972b]. The analytic 
solution of this paper carried out for a spherical Moon shows 
that the front side spherically symmetric solution is more accurate 
than the results of Reisz et al . [1972] would indicate. 

The solar wind distorts the dipole and higher order multi- 
pole parts of the time dependent induced lunar field in a manner 
somewhat similar to the distortion of the Earth's dipole field 
by the solar wind. Thus we speak of a time dependent lunar 
magnetosphere; its boundaries are delineated by the front side 
lunar surface, the diamagnetic cavity walls and the interior 
lunar bulk electrical conductivity profile [Sonett and Colburn, 
1968; Johnson and Midgley, 1968]. The sunward surface of the 
Moon corresponds to the Earth's magnetopause. In the case of 
the Earth's field, the magnetopause is formed where the magnetic 
pressure of the dipole field is comparable to the dynamic pressure 
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of the incoming plasma. In the lunar case, the magnetic 
pressure of the induced fields is insufficient to balance the 
solar wind dynamic pressure. Rather, the lunar surface, which 
is a fixed boundary and absorbs much of the incident plasma, 
marks the sunward extent of the induced field. The lunar 
magnetosphere concept is supported by the theoretical results 
presented here and the empirical data of Part II. For example, 
in this paper we show that for a lunar magnetic field induced 
by an interplanetary field transverse to the direction of solar 
wind flow, there are two points on the surface of the sunlit 
hemisphere where the magnetic field tangential to the surface 
vanishes. These points are qualitatively similar to the neutral 
points found on the Earth's sunward magnetopause [Spreiter and 
Alksne, 1970]; however, for the Moon the positions of vanishing 
tangential field are not neutral points because the radial field 
does not vanish there. 

The next section gives a description of the theoretical 
model. Following that we proceed to a discussion of the re- 
sults of the analysis. Mathematical details are given in the 
appendices. This will permit the reader mainly interested in 
the configuration of the distorted induced lunar magnetic field, 
i .e. the lunar magnetosphere, to avoid the mathemati cal complexity 
of the theory. 
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THEORETICAL MODEL 


Induced lunar magnetic field lines must accommodate them- 
selves to an available volume delineated by the front side 
lunar surface, the cavity boundary and> the lunar electrical 
skin depth profile. A precise radial, location of the internal 
boundary to this volume cannot be determined in general since 
attenuation of electromagnetic waves occurs over a range of 
depths. However, an equivalent depth can be defined for a 
given frequency, such that a correct measure of the inductive 
response is obtained by considering the Moon to be nondissipa- 
tive above this depth and perfectly conducting below. Thus the 
exclusion of fluctuating fields from the lunar interior as a 
result of cumulative attenuation is modelled with a perfectly 
excluding core of an equivalent size. The model adopted here 
is thus a quasistatic one in which the induced lunar field lines 
are distorted to occupy the space between an infinitely con- 
ducting core of radius b and the confining boundaries of the 
front side lunar surface of radius a and the walls of the cylindri 
cal cavity of radius a which extends downstream from the Moon 
(positive z-di rection) . Both the lunar shell of thickness (a-b) 
and the downstream cavity are nonconducting. 

The center of the Moon coincides with the origin of the 
coordinate system and the axis of the cylinder is the z-axis. 

We use both spherical coordinates (r, 9, cp) and cylindrical 



coordinates (p, cp, z). The geometry is illustrated in Figure 1. 
The space outside the Moon and the cavity is occupied by solar 
wind or magnetos heath plasma, and magnetic fields, associated 
with induced currents in the Moon, on the lunar front side 
surface, or on the cavity boundary, cannot penetrate this 
region. 

At the frequencies of interest, the lunar cavity radius 
is much smaller than the vacuum wavelengths of electromagnetic 
waves, thus displacement currents are negligible. In the non- 
conducting shell of the Moon, and in the cavity, the magnetic 
field is derivable from a scalar potential Q 


H = vn , ( 1 ) 

7 2 Q = 0 . (2) 

The theory assumes that the external fluctuating field 
is spatially uniform. Most of the frequencies of interest 
are sufficiently low so that the solar wind sweeps through an 
appreciable distance in times comparable to the periods 
associated with the induction. For example at a frequency 
of 10 Hz, a relatively high frequency for lunar surface 
magnetometer observations, and a solar wind speed of 400 km/sec, 
the solar wind travels 4x1 0^ km or about 20 lunar radii in a 
wave period. Because of the cavity deformation of the induced 
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field, dipole and higher order multipoles of the induced field 
are present; however, the high frequency modification of these 
multipoles, due to spatial nonuniformity in the external field 
which is present at sufficiently short time scales, is not in- 
cluded in the theory. 

Everywhere external to the Moon and cavity we assume the magnetic 
field to be uniform, of magnitude unity (without loss of 
generality), and oriented either normal to the cylindrical axis 
e.g. in the y-direction, or parallel to the cavity axis, i.e. 
in the z-direction. By superposition, the solution for an 
arbitrarily oriented external magnetic field can be obtained 
from the fundamental solutions presented here. The boundary 
conditions to be satisfied by the solutions of (2) are conti nui ty 
of the normal component of magnetic field on the cylindrical, 
boundary and the front side lunar surface and zero magnetic 
field normal to the surface of the inner core. In addition, 
the magnetic field in the cavity and in the insulating shell 
must be continuous on the far side lunar surface. 

The possible effects of diamagnetism on the induced 
field are not included. Diamagnetism influences the steady 
field by leading to an increase of the steady interplanetary 
field inside the cavity. We are concerned here only with the 
time varying field, but it also might be influenced by fluctua- 
tions of the intensity of diamagnetism as the solar wind para- 
meters vary. We also neglect the effects, other than confinement. 
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of the undoubtedly complex processes occuring in the interaction 
layer of the cavity walls. 



DISCUSSION OF RESULTS 


The details of the analytic solution are presented in the 
appendices. Appendix A deals with. the case of the external 
magnetic field perpendicular to the cavity axis and appendix 
B considers the external field aligned along the cavity axis. 

For each orientation of the external field we write the scalar 
magnetic potential and the fields in the nonconducting lunar 
shell in terms of expansions in spherical harmonics. Similarly, 
in the cavity, the. potential s and fields are given by expansions 
in terms of cylindrical eigenfunctions of the Laplace equation. 
This facilitates the application of the boundary conditions on 
the cavity walls. The solutions in the cavity are then written 
in terms of spherical harmonics with the aid of an identity 
which permits expansion of cylindrical eigenfunctions in terms 
of spherical ones. Finally, application of the appropriate 
boundary conditions on the lunar surface r=a determines the solu- 
tions . 

In addition to the asymmetric solutions, formulas for the 
magnetic fields in the spherically symmetric vacuum SSV approxi- 
mation (A. 10) - (A. 13) and (A. 48) - (A. 50) and the spherically 
symmetric plasma SSP approximation (A. 6) - (A. 9) and (A. 45) - 
(A. 47) are given in the appendices. In the SSV case, the field 
is just the superposition of the uniform external field and an 
induced dipole which prevents the external field from penetrating 
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the infinitely conducting core. In the SSP case, the solution 
is similar except for an amplification of the strength of the 
induced dipole as a result of compression of the induced field 
into the lunar shell. 

If the radius of the lunar core is zero the magnetic field 
is everywhere the uniform field existing outside the space of the 
Moon and cavity. This may be seen in the case of the external 
field normal to the cavity axis, for example, from (A. 28) and 
(A. 29). With b=0 , a ^ =1 and the equations for the coefficients 
are homogeneous. Thus C^=0 and from (A. 15) " (A. 17) it 
follows that the magnetic field in the cavity is simply the 
uniform external field. Further, from (A. 30) and (A. 31) it is 
clear that Aj-1 and A^=0 for 1*2. Using (A. 5) we find B^=0 
for all t . Thus the field inside the Moon is indeed the 
uniform field of the external space. 

In the following we present results showing magnetic 
fields and transfer functions for different external field 
orientations and core radii between 0.5 and 0.9 lunar radii. 
Differences between the spherically symmetric approximations 
and the exact solutions of this paper generally become increas- 
ingly significant as the core radius increases, corresponding 
to waves whose attenuation is progressively more significant at 
shallower depths in the Moon. Thus the inadequacy of spherically 
symmetric approximations will be most serious at the high frequency 
end of the spectrum accessible with the lunar surface magnetometers. 
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Lunar Surface Fields for External Field Parallel to Cavity Axis 

Figure 2 shows the radial component of the total magnetic 
field on the dark side lunar surface (r=a), O<0< tt/2 , as a 
function of 9 . Also shown is the radial magnetic field computed 
from the SSV approximation. All fields are scaled to an ex- 
ternal field of unit amplitude. On the sunward side of the Moon 
the radial magnetic field is identical to the radial component 
of the external field, cos0 . In the SSV approximation the 
radial field increases monotoni cal 1 y from zero at the limb to a 
maximum at the antisolar point and shows a cos0 dependence. 

Also, as the core radius increases, for a given 9 , the radial 
SSV field decreases. For all b/a the radial field is larger 
than that of the SSV approximation, the difference becoming 
more significant as b/a increases. For b/a = 0.5, 0.6, and 0.7 
the radial field increases monotoni cal 1 y from zero at the limb 
to a maximum at the antisolar point. However, for b/a = 0.8 
and 0.9 the radial field increases to a maximum very near to the 
limb and subsequently goes through a minimum as 0 approaches 0°. 
As a result, near the limb the field decreases with increasing 
core size, whereas near the antisolar point the field increases 
with increasing core size. The increase in the radial field 
above the SSV value, the local minimum in the field seen near the 
antisolar point and the field maximum near the limb for large 
cores can all be understood in terms of the distortion of the 
induced field which adjusts to the available space defined by the 
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nonconducting shell and the cavity. For the larger values of b/a 
the combined fields of the cavity and interior currents add to 
raise the total to a very large value near the terminator. Away 
from the terminator, i .e. towards the antisolar point, the total 
field drops quickly but never to a value equal to the SSV case. 
Thus either the cavity edge field is of very high order and 
therefore of short range, or its direction is rapidly changing 
so that its contribution to the radial component of the lunar 
far side surface field attenuates rapidy as 0-*O, or both these 
effects are present. 

Figure 3 shows the 0-component of the total field at the 
surface of the Moon normalized as before to unity driving field 
and the 0-component of the SSV field for varying b/a. The SSP 
fields on the front side lunar surface are virtually identical 
to the fields of the exact solution. Front side amplification 
of this field is thus negligibly affected by the day-night 
asymmetry of the lunar plasma environment. On the lunar far 
side, Hq is also amplified considerably above its SSV value, 
amplification increasing with increasing b/a and proximity to 
the limb. As one might expect the effect of cavity confinement 
of the far side field dominates near the terminator where Hq 
attains a maximum. 

It is interesting to note that the calculations show the 
maximum to move slightly forward of the terminator with increas- 
ing b/a. For b/a = 0.9 this maximum is about 3-4 degrees ahead 
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of the limb (a precise location of the maximum would require 
evaluating the field on a finer 0 scale than was done for this 
paper). This is a consequence of field line redistribution 
internally in the Moon as a result of the day-night asymmetry 
of the space external to the Moon. Figure 3 clearly shows a day- 
night asymmetry in Hg which becomes more pronounced for increas- 
ing b/a. This is readily understandable since the compressive 
effects of the solar wind are different from those of confinement 
by the cavity. Lastly we note that for the case of the external 
field aligned along the cavity axis vanishes and H r and H Q 
are independent of cp because of azimuthal symmetry. It is 
remarkable that for the external field aligned case the front 
side lunar surface magnetic field is practically identical to the 
SSP approximation while the dark side surface field is significantly 
different from the SSV field even for large b/a 1 

Lunar Surface Fields for External Field Perpendicular to Cavity Axis 

We next turn to the configuration where the external magnetic 
field lies normal to the cavity axis. Figure 4 shows the radial 
component of the field on r=a vs. 0 between the antisolar point 
and the terminator together with the SSV approximation to this 
field. At the limb H r must be unity to preserve continuity of 
the normal component of the magnetic field. Similarly on the 
sunward side of the Moon this continuity condition prescribes 
the variation of H r with 0 to be that of the external field. We 
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note that H r «= sincp so that the variation of H r shown is for the 
y-z plane. Alternatively, for an external field in the x-direction 
the variation of H r in Figure 4 would be for the x-z plane. 

The dependence of H r on 0 is very different from that of 
the SSV approximation. Near the antisolar point the asymmetric 
theory gives values of H r which are smaller (increased induction) 
than those found from the SSV theory. However, the situation 
reverses as the terminator is approached where the values of H- 
are larger than the SSV values. The changeover takes place at 
a value of 0 determined by the size of the core; as b/a is in- 
creased the crossover occurs nearer the limb. Further the SSV 
radial functions are concave downwards everywhere on the dark 
side while the asymmetric calculations show radial functions 
which display upward concavity over large 0 intervals. It is 
particularly interesting to note the case of b/a = 0.9 where H r 
for asymmetric induction is nearly zero over a range of 0 up to 
60° from the antisolar point. This can be understood by noting 
that the distortion of field lines around the infinitely con- 
ducting core must be accommodated entirely within the cavity 
space, while in the SSV approximation the field line distortion 
is more gradual and can occur in the entire space external to 
the core. 

Figure 5 shows the calculations for Hg which is « sincp (see 
Appendix A). This figure views Hg in the plane cp = tt/2, i.e. 
the plane defined by the external field and the cylinder axis. 
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For convenience we can think of the forcing field to lie in the 
plane of the ecliptic, which then means that Hg is measured 
equatorial ly about the Moon. It is clear however that the plane 
as defined above can take on any orientation that includes the 
cylinder axis. For comparison, Figure 5 also shows the SSV dark 
side field for b/a = 0.9 and the SSP front side fields for b/a = 
0.8 and 0.9. The fields on the front side are generally smaller 
than the SSP fields. For b/a = 0.9 the differences between 

actual and SSP fields are largest and amount to nearly 20 per 
cent at the subsolar point. The differences diminish as b/a 
becomes smaller. On the dark side Hg is larger than the SSV 
approximation. There is considerable cavity amplification of Hg 
which can increase as the limb is approached from the antisolar 
point. The SSV approximation always yields a monotonic decrease 
of Hg with proximity to the limb. 

An important aspect of the distortion introduced by the 
asymmetry is that H 0 vanishes for a value of 9 > 90°, the zero 
of Hg occuring more forward of the limb for increasing b/a. This 
is qualitatively similar to the case of the Earth's magnetosphere 
as discussed in the introduction. Generally speaking the dynamic 
pressure of the solar wind is so much greater than the pressure 
of the induced field that one should think in terms of field 
line redistribution here rather than sweeping back of field lines 
to the dark side by the action of the solar wind, as is often 
mentioned as the mechanism for the formation of the Earth's 



- 16 - 


magnetic tall. The case b/a = 0.9 shows a nonmonotonic varia- 
tion in Hg on the dark side of the Moon which is attributable 
to the 8-component of the summed fields of the cavity and core 
currents peaking at an angle 9<90°. This effect essentially 
disappears for small b/a. The decrease in the sunward response 
of the Moon from that of the SSP approximation is attributed to 
redistribution of field lines from the front to the back. 

A discussion of the field is not presented here since 
it will be seen that can be viewed as being identical to the 
North-South transfer function which is discussed in a later 
section. 

Magnetic Fields on the Cavity Boundary 

We have computed the magnetic fields along the walls of the 
idealized cavity for the two cases of the external field along 
and normal to the cavity axis. In the next two figures we use 
cylindrical coordinates to describe the magnetic field. Figure 6 
shows the field H z on the cavity wall for the case of aligned 
external field as a function of distance along the cavity from 
the terminator. The abscissa also shows the angle 8 correspond- 
ing to the varying distance. Since the normal component of 
magnetic field is required to be continuous on the cavity wall, 

= 0. Also because of azimuthal symmetry vanishes. The 
field is a maximum at the terminator and decreases with distance 
along the cavity, becoming vanishingly small for distances in 
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excess of about one lunar radius downstream from the limb, 
corresponding to 8 of « 45°. Also H z increases in magnitude 

for larger b/a . Azimuthal surface currents of density K = H 

cp Z 

(MKS units are understood) flow on the cavity wall and are 
responsible for the cavity confinement of the induced fields on 
the dark side. As in the case of the field, the current density 
attains its maximum value at the terminator. In the neighbor- 
hood of the limb the cavity wall currents are of the same order 
of magnitude as those which flow on the sunward side in response 
to the confinement there of the induced magnetic field. The 
surface current density on the sunlit hemisphere can be found 

from K = - s i n9 - H Q . 

cp u 

We turn next to the case of the external magnetic field 
perpendicular to the cavity axis. Figure 7 shows the magnetic 
fields H z and along the cavity walls for this case. The 
field component H p is simply that component of the external 
field by continuity of H p at p=a. Note that « coscp and 
H « sincp. Thus we view H in the plane which includes the 
cavity axis and the normal to the external field, and we view 
H z in the plane which includes the cavity axis and the external 
field. The field H z peaks at about 0.5 lunar radius downstream 
from the terminator. The importance of the ratio b/a in determin- 
ing the magnitude of the cavity wall fields is apparent. In 
this case two current components coexist on the cavity walls. 
Azimuthal currents of density = H z are found as in the aligned 



-18- 


field case and we find this circulating current system peaking 
downstream at 0.5 lunar radius. Additional currents flow in 
the z-direction with density K z = - coscp. The current system 

on the cavity boundary is of course continuous at the terminator 
with the current system on the sunward side of the Moon which is 
given by = cos0 sincp - H 0 and Kg = - + coscp. No dissipa- 

tion arises from these currents and they are to be regarded as 
Hall currents flowing across the local field. Currents flowing 
across the front surface of the Moon to conform with the con- 
finement requirement close in a complex pattern over and under 
the lunar cavity. Cavity compression of induced fields may be 
attributed to the currents K , K on the cavity boundary. 

2 Cp 

Transfer Functions 

For comparison to experimental data as will be done in 
Part II it is convenient to describe the lunar response in terms 
of transfer functions, i.e. ratios of the observed field com- 
ponents to corresponding components of the external field. The 
empirical data presented in Part II is given in terms of radial 
or local vertical, North-South and East-West transfer functions. 
The directions are all local at the site of the Apollo 12 surface 
magnetometer wh i ch is essentially in the ecliptic plane. Con- 
tributions to the North-South transfer function can only be made 
by an external field which is perpendicular to the ecliptic 
plane. Thus the North-South transfer function is identical to 
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the field viewed in the ecliptic plane and driven by an 
external field of unit amplitude which is normal to the cavity 
axis and the ecliptic plane. The North-South transfer function 
is independent of the absolute spatial orientation of the ex- 
ternal field as long as it has some nonzero component normal 
to the ecliptic plane. 

Figure 8 shows the North-South transfer function calculated 
from the asymmetric theory as a function of 9 together with the 
SSV and SSP values for comparison. Both the SSV and SSP approxi- 
mations to this transfer function are independent of solar 
longitude (solar longitude or local time at the site of the 
surface magnetometer is identical to 0 in this view of the 
North-South transfer function). For the SSP approximation we 
find from (A. 8) or (A. 9) that the tangential transfer function 
i s 


and for the SSV approximation it follows from (A. 12) or (A. 13) 
that the tangential transfer function is 



The effects of asymmetric confinement stand out clearly. As in 
the discussion of the fields, the larger b/a, the greater the 
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departure of the asymmetric transfer function from the spherically 
symmetric approximations. We note that confinement is strongest 
on the sunward side of the Moon, though not as great as that pre- 
dicted by the SSP approximation. The reduced amplification on 
the front side below the SSP value can be understood in terms of 
the leakage of induced field lines into the cavity space and the 
redistribution of front side field lines to the back with an 
accompanying decrease in front side compression. Amplifications 
are in agreement with the values found experimentally (see Part II). 
The values of the transfer function decrease monotoni cal 1 y as one 
moves from the subsolar towards the antisolar point indicating 
that effective confinement becomes less perfect as the dark side 
of the Moon is approached. Nevertheless, it can be seen that in 
spite of the assumption of perfect vacuum in the lunar cavity, 
there is substantial amplification of the North-South magnetic 
field on the dark side, where values of amplification are always 
larger than those of the SSV approximation. This increase in the 
far side amplification above the SSV value can be understood in 
terms of cavity compression of the induced field. For b/a > 0.8 
the amplification is everywhere greater than the theoretical high 
frequency limit (1.5) for a conducting Moon in vacuum. 

Unlike the North-South transfer function, the radial and 
East-West transfer functions depend on the absolute orientation 
of the external field. External fields in the ecliptic plane 
induce magnetic fields which determine these transfer functions. 
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Since these external fields can be arbitrarily oriented in the 
ecliptic plane, the radial and East-West transfer functions 
cannot be a-priori defined. Their values depend on the relative 
contributions of external field components parallel and perpendicu 
lar to the cavity axis. 

As an example we show the radial transfer function for an 
external field which is at 45° to the cavity axis and in the 
ecliptic plane and co = 90°. Figure 9 shows this transfer function 
for varying b/a as a function of 9. The boundary condition on 
the normal component of the magnetic field requires that the 
radial field component be continuous across the front surface of 
the Moon so that on the entire front side this transfer function 
is unity for any b/a. On the dark side of the Moon a minimum is 
seen at about 30 degrees away from the antisolar point. As 
before the variation in the transfer function increases with b/a. 
Note that the transfer function in the SSV approximation is in- 
dependent of solar longitude (as in the case of the North-South 

transfer function solar longitude or local time is identical to 

b 3 

8) and is given by 1 - — ^ (see (A. 49) ). 

a - 5 

At the subsolar and antisolar points, Hg , H and the 
tangential (North-South and East-West) transfer functions are 
identical . Only external fields normal to the cavity axis 
induce fields which contribute to these quantities. Figure 10 
shows the tangential transfer functions at the antisolar and 
subsolar points together with the appropriate spherically 
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symmetric approximations vs. b/a . At the subsolar point the 
transfer function is smaller than the SSP approximation, the 
difference being negligible for small cores, i .e. low frequency 
waves, and important for large cores, i.e. high frequency waves 
which damp at relatively shallow depths. At the antisolar point 
the transfer function is larger than the SSV approximation, the 
difference becoming more significant as the core size increases. 

The variation with b/a of the radial transfer function at 
the antisolar point is presented in Figure 11. The SSV approxi- 
mation to the transfer function is consistently smaller than 
the radial transfer function of the asymmetric theory. For a 
given lunar conductivity model the SSV approximation over- 
estimates the induction. 

Finally in Figure 12 we show the North-South transfer 
function at the limb vs. b/a. The radial transfer function at 
the limb is unity and the East-West transfer function at the limb 
is dependent on external field orientation. Neither the SSP 
nor^SSV approximations to the North-South transfer function are 
a-priori expected to be good at the terminator. Both approxi- 
mations are also shown in the figure. The transfer function has 
values intermediate between the SSV and SSP approximations. The 
North-South magnetic field shows considerable amplification at 
the limb with values closer to those of the SSP than to those 
of the SSV approximation. 
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CONCLUSIONS 

The analytic solution of this paper predicts the ways in 
which the confinement attributable to the diamagnetic cavity 
influences lunar magnetic induction and permits a quantitative 
assessment of the importance of these effects for various 
frequencies and for both sunward and dark hemisphere transfer 
functions. The results will be used in Part II to compare the 
theory with experimental transfer functions for both sunward 
side and dark side data. 

Field line leakage into the cavity and the internal re- 
distribution of field lines from the front side of the Moon 
to the night side contribute to a reduction of the sunward side 
field line compression from that which would be predicted by 
the spherically symmetric plasma approximation. Cavity com- 
pression of field lines, consistent with surface currents on the 
cavity boundary, acts to amplify the tangential magnetic fields 
and significantly alters the radial magnetic field on the lunar 
dark side from the prediction of the spherically symmetric 
vacuum approximation. These effects are most important for waves 
whose equivalent exclusion depths are small. 
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F i gure 1 . 


Figure 2. 


FIGURE CAPTIONS 

Geometry of the theoretical model. The Moon has 
radius a with a spherical core of infinite electri- 
cal conductivity, , surrounded by a nonconduct- 

ing shell of thickness (a-b). A cylindrical non- 
conducting cavity of radius a extends downstream 
from the Moon. The coordinates r, 9 , cp are spherical, 
p, cp, z are cylindrical coordinates, with the origin 
of both coordinate systems at the center of the Moon 
and the z-axis coinciding with the cylinder axis. 

The external field is uniform and is either normal 
to the cavity axis, i.e. in the y-di rection (cp = tt/2) 
say, or parallel to the cavity axis. The variable 
core radius allows investigation of waves of different 
frequency which are excluded from different portions 
of the Moon. 

The magnetic field component H r and the SSV approxi- 
mation to H r as a function of position on the dark 
side lunar surface, r=a , O<0< 90°. The external 

' . - , / • ’ : tv 

field has unit ampl i tude and is para,l l, el., to the 

* ’ > * 

cylinder axis. The antisolar point is 9=0° while 
9=90° is the limb. Note the increase in the radial 
field above the SSV value, the local minimum in the 
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Figure 2 (continued) 

field near the antisolar point and the field maximum 
near the limb, the latter two phenomena for large b/a . 

Figure 3. The magnetic field H 0 on the lunar surface as a 

function of 0 for various b/a. The external field 
has unit amplitude and is aligned with the cavity 
axis. The field according to the SSV approximation 
is shown for b/a = 0.5 and 0.9. The field from the 
SSP approximation is virtually identical to the asym- 
metric field for 90°<9<180°. There is considerable 
amplification of H 0 on the lunar dark side. Further 
the amplification is asymmetric with respect to the 
terminator. Both the magnitude of the amplification 
and the extent of the asymmetry in the lunar response 
increase with b/a. 

Figure 4. The radial magnetic field on the lunar dark side 
surface with b/a as a parameter when the external 
field of unit amplitude is aligned normal to the 
cavity axis. The radial field according to the SSV 
approximation is also shown. On the lunar front side, 
H r is identical to that component of the external 
field. The radial field is viewed in the plane of 
the external field and the cylinder axis. The field 
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Figure k (continued) 

is reduced below the SSV value near the antisolar 
point but exceeds the SSV value near the terminator. 
For b/a = 0.9 the field is essentially zero within 
60° of the antisolar point. 

Figure 5. The magnetic field component Hg on the lunar surface 
vs. 0 for various b/a in the case of unit external 
field normal to the cavity axis. The field according 
to the SSV approximation for b/a = 0.9 is 

shown for 0<9<90° and Hg from the SSP approximation 
for b/a = 0.8 and 0.9 is shown for 9O°<0<18O°. The 
field is viewed in the plane of the external field 
and the cavity axis. On the front side near the sub- 
solar point Hg is reduced below its SSP value, while 
on the dark side near the antisolar point Hg is 
amplified above its SSV value. The field has a zero 
forward of the terminator as a result of the redistri- 
bution in the lunar interior of field lines from the 
sunward to the dark side hemisphere. 

Figure 6. The field H z on the cavity boundary p=a as a function 
of distance downstream from the limb with b/a as a 
parameter. The unit external field is parallel to 
the cavity axis. Note the attenuation in H z with 
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Figure 7. 


Figure 8. 


distance from the limb; beyond about one lunar 
radius H z is essentially zero. The figure also 
shows the distribution of azimuthal surface current 

dens i ty K = H . 

' cp z 

The field components H z and , or alternatively the 

surface current densities K = H and K = H - 1, 

cp z z cp 

on the cavity boundary p=a as a function of z for 
various b/a. The unit external field is perpendicular 
to the cavity axis. The currents K or fields H 

' Cp 2 

decay monoton i ca 1 1 y with distance downstream from 

the terminator. The currents K or fields H attain 

z cp 

a local maximum downstream from the limb before decay- 
ing to vanishingly small values. 

The North-South tangential transfer function vs. solar 
longitude or local time with b/a as a parameter. The 

curves also represent H for a unit external field 

r cp 

perpendicular to the cavity axis if H is viewed in 
the plane containing the cylinder axis and the normal 
to the external field. Values for the SSP and SSV 
approximations are also shown. For clarity of the 
figure, the b/a parameter has been omitted from the 
SSV lines. These values are b/a = 0.6, 0.7, 0.8 and 
0.9, proceeding in order from the lowest to the highest 
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Figure 8 (continued) 

SSV line. The transfer function or H is inde- 

cp 

pendent of 0 in the spherically symmetric approxi- 
mations. Note the reduction in front side amplifi- 
cation below the SSP value and the increase in the 
dark side amplification above the SSV value. 

Figure 9. The radial transfer function vs. solar longitude for 
several b/a. The external field is at 45° to the 
cavity axis and in the ecliptic plane and cp = 90° . 

For the entire front side, the radial transfer function 
is unity. On the lunar dark side, the transfer 
function shows a minimum about 30° from the antisolar 
point. The transfer function in the SSV approximation 
is independent of solar longitude. 

Figure 10. The fields Hq , or alternatively the tangential 

transfer function at the subsolar point 0 = 180° and 
at the anti solar point 0=0° as a function of core 
size. The cavity amplification of fields at the anti- 
solar point and the reduced amplification of fields 
at the subsolar point due to field line leakage into 
the cavity become increasingly important for larger 
cores, corresponding to waves whose effective exclu- 
sion depths are progressively shallower. 
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F i gure 


Figure 


11. The radial transfer function at the antisolar 
point vs. b/a together with the SSV approximation. 

For a given lunar conductivity model the SSV approxi- 
mation overestimates the induction. 

12. The variation with b/a of the North-South transfer 
function at the limb. Both SSV and SSP approxi- 
mations are also shown. The North-South magnetic 
field shows considerable amplification at the limb 
with values closer to those of the SSP than to those 
of the SSV approximation. 
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APPENDIX A 

ANALYTIC SOLUTION FOR EXTERNAL MAGNETIC 
FIELD NORMAL TO THE CAVITY AXIS 

The potential in the insulating shell Q* is 

00 

n 1 = Y, ( A £ pt + B -i r "^ _1 ) P 1 ( cos0 ) sinc P » (A.i) 

t-i 

where A. and B. are coefficients to be determined as part 

-V \/ 

of the solution, and p|(cos9) are the associated Legendre 
functions of order 1 (the sign convention used here is such that 
for example pj(cos9)= sin9). The cp dependence of O 1 follows 
from that of the external magnetic field which is 


A A ^ 

r sin9 sincp + 9_ cos9 sincp + £ coscp , 

or 

A A 

p s i ncp + cp coscp * 


AAA A A A 

where r_, 9., £ and _p, £, z are the unit vectors of the 
spherical and cylindrical coordinate systems, respectively. 
The magnetic field in the shell is 

oo 

H* = Y I'LA^r^'” 1 -(£+1 )B^r"^" 2 } p|( cos9 ) s i ncp 
1 = 1 


(A. 2) 
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H e ■ Z 

1 = 1 


f A ^-1 , D 2) dP l (co s 0 ) . 
( A i r +B i ) — are s,ncp 


(A. 3) 


u i _ r / A 

cp U. \ A l 

1=1 


.t-1 , „ p l< cose > 

• d o r / ■ ■ t' a 


/ s i ne 


coscp 


(A. 4) 


The coefficients A^ and are related by the condition 
that the radial component of magnetic field be zero on r = b 


l A^' 1 - U+l ) B^b"' 1 " 2 = 0 


(A. 5) 


The normal component of magnetic field must be sin9 sincp on 
r=a,^^0^n . For the case of spherically symmetric 
plasma conf i nement (SSP) , this condition is applied over the 
entire range of 9, i .e. 0 to tt . The solution is then straight" 
forward with 


A 1 " 


,3 k 3 


a -b 


1 — 7 » B l = ~' q , b tt > A, , B p = 0 for t z 2 
u a J -b J 


'l* l 


(A. 6) 


and in the insulating shell, region I , 


H = a ^ s '^ 9 sincp f ^ b^ N . 

r ( a ^- b 3 ) V ' T 7 


(A. 7) 


u _ a i cos 6 sincp / , b* \ 

He ■ (.i. b j) ° + i?> 


(A. 8) 
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H * = 


|22S- (l + 

(a 3 -b 3 ) V 2r 3 ' 


In the spherically symmetric vacuum approxi mat ion (SSV) i t 
is simple to verify that 

3 

Aj = 1 , Bj = , A^, = 0 for t s 2 , 

and in region I 


f-L = 


• - • / « D \ 

= s i n0 s i ncp (J - —jj 


H 0 = 


- 


, b 3 , 

co s 0 si ncp (1 + « j 

v 2r^' 


coscp (j + ^-j) 
2 r -' 


The potential in region II, the cavity, is given by 


00 , 

IT ~ “k z 

Q 1 = psincp + si ncp > C e n J,(k p) 

i—i n in 

n=l 


where C p are coefficients to be determined in the solution, 

J, is the Bessel function of order 1 and k a is the nth zero 
I n 

of jj ( k p a ) = 0 (the prime denoting the derivative). The 
magnet i c field is 



s i ncp + s i ncp / 

n=l 


C n e 


■k z 
n 


9 


(A.9) 

(A. 10) 

(A. 11) 
(A. 12) 
(A. 13) 

(A. 14) 


k p J((k n o) 


(A. 15) 
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H** = coscp + coscp ^ C p e n p ^ J^(k n p) , ( A . 1 6 ) 

-<P n=l 

oo |< Z 

H** = - sincp ' C n k n e ° * (A. 17) 

2 n=1 

This sol ution automatical 1 y satisfies the continuity of H p 
on the cavity boundary p = a since jj (k n a) = 0* To apply 
the boundary conditions on the hemispherical surface r = a , 

0 <; 0 < rr/2 , (A. 14) - ( A - 1 7 ) are most conveniently rewritten 

in terms of spherical coordinates. With the aid of an identity 
due to Cooke [1956] we find 


“ “ ( k r) n ( - 1 ) n+ ^ 

Q 11 = rsine sincp + /. / — pI( cos0 )sincp, (A. 18) 


£l * n=l < n+1 > i 


H 


II 


, n+1 


Sine sincp + V C, 4^ 1 ^ -( k » r )" p'(cos0),(A.19) 

“ 4 r n ti (n+1)!' 1 " 


t = ] 


H 


II 


r-> t i \ n +l n d P "L ( COS 9 ) 

cose sincp + £ C t ^ k i r) 


4 -l * r n=l 


(A. 20) 


H 


II 

■cp 


V' rn # 1 \I"I + 1 n P n (cOS0) 

\ „ coscp V (-1) /i. ,\ii n 

coscp + ^ Z t -p— l ^ ( nil )'! (k £ r) “TTrTS 

1=1 n = l 


(A. 21) 
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Continuity of H r on r = a yields 

00 

£ p|(cos9) - sin9 = 

4 = 1 

1 r 1 i <iiC? < k t a >" P n < cos0 > > 0<e " /2 

4 = 1 n=l 

or 0 , y <0< n 

Equation (A. 22) leads to the following relationship among the 
unknown coefficients 


( 24 '+1 ) L 1 1 J 


4 

3 


V ,1 


= ) — — Y LlU n fk a) n P / l'=] 2 3*** 

A “ L , M \ []< l a) P n4 ’ ^ 1 ’ z ’ 3 

4-1 n = l (n+1) 1 


where 6^, , ^ - 1 if 4 ' = 1 

=0 if l' 4 1 


P 44' 


p[ 2 (0)^ 


(0) 


1 dp !' 
p{ (0) ^ 


( 0) 


; r* 4.1 ^ 


4 , (4 / +l ) 
(24'+l ) 


, 1 = 1 ' , 


(A. 22) 


(A. 23) 


(A- 24) 


(A. 25) 
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In (A. 25) the argument of the associated Legendre functions 
and their derivatives is cosQ = 0 . 

On r=a, O<0< tt/2, H q and are continuous. It is 
sufficient to apply one of these conditions, the other being 

redundant. Continuity of H leads to 

cp 


V 

/ 

i—x 

t = 1 


(A^a^ -1 + B^a"^" 2 ) p|(cos0) = sin0 

(*>" P^tcose, 


(A. 26) 


for O<0< tt/2 . To generate a set of equations for the 
unknown coefficients from (A. 26) it must be noted that on the 
range of 0 given the odd associated Legendre functions form 
a complete set. Expanding (A. 26) in terms of this set of odd 
associated Legnedre functions we find 


l k 

t = 1 


t - 1 




-i-v 


p n' 


2 

I 


V , i 


oo 



1=] 



a 


V' 

L 

n = l 


i i \ n+ 1 

(nil’)'!' (k t a) P nl 


(A. 27) 


for = 1 ,3,5,*** 

Equations (A. 5), (A. 23) and (A. 27) suffice to determine 
A^ ,B^ and . After considerable algebraic manipulation, a 
single set of equations for may be derived 



8 I/OH 
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V 


... 1 


c r (k a) v i'u'+') (4- - n 

(_-, {_a 


(2t'+l )(t'+ 1 )! p! /(0) 


dP 


I 


n even 


( k m a ) n < 0) ( ? (o V + a n ) ' 1 } 

( n+1 ) I l'(t'+\ )-n(n+l ) 


L 

n odd 


(k a ) n n pJfO) r 
v m ' n v J_ \ 


(n+1) 1 


t even 


a^( 21 +] ) 
21 ( 1 +] ) 


,dp! x2 

(0) ) 


[l(l+ ])-l’(l’ +] )][-t(t+l )-n( n+1 )_ 


-} 


= -j ( i - a i ) V j 


l' = 1 ,3,5, 


where 


21 +] 

1 1 ITT (|) 

a t " , . 2-t+K 

<1 - (|) ) 


Once is determined from (A. 28), A^ may be found from 


(A. 28) 


(A. 29) 
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V 


1 -] 


l a® < k m a > 

m=l 

TT+T 


i l ’ ] 

+ 577 


2 (,- ( |, )(t +1); '-(f) 


(2t+l )p!(0) 


V m V 

o / k 2T^fT a _Z-» 

2t 2 (*.+l)(l-(|) ) m "' 


n even 


dP n 

n(k m a) (0) 

( n+1 ) i “ 1(1+1 )-n( n+1 )| 


for l odd and 


}' ^(ka ) 1 

L a v m ' 


V 


t - 1 m=l 


2(1-(!) M + ')U+1) 


dp! 

( 2 ^’> gjsF <°> 

+ - - 2-t + l v 

2l Z (l + ]){]- (§) ) 


V C 

y ^ 

itt! a 


y 


m 


n odd 


n(k a) n pI( 0) 
m n 


(n+1 )! ) - n( n+1 )} 


for l even. With and determined, follows from (7). 
The zeros of Jj are tabulated in Abromowitz and Stegun [1964] 


(A. 30) 


(A. 31) 
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The infinite set of algebraic equations has been truncated 
and solved numerically. Convergence of the solution has been 
assured by increasing the size of the set until a further 
increase in size produced a negligible change in the results. 
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APPENDI X B 


ANALYTIC SOLUTION FOR EXTERNAL MAGNETIC 
FIELD PARALLEL TO CAVITY AXIS 


The external magnetic field is 


AAA 

z = cos9 r - sin0 0 


Thus the magnetic potential in the lunar shell is 


O' 


1 (v* + B i r ‘ t ‘ 1 ) P^( cose ) 


(A. 32) 


1 = ] 


where P^Ccos©) are the Legendre polynomials and the magnetic 
field ' s g i ven by 


H r * I («t 

t = 1 


r l ' ] - (1 + ] )B, r' 1 ' 2 ) P^(cos0) 


(A. 33) 


H 9 " ), (a^ 1 + B^r 

l=\ 


-l-2\ dF V cos0 ) 




(A. 34) 


With the azimuthal symmetry characteristic of this external 
field alignment, all quantities are independent of cp and there 
is no azimuthal magnetic field. The coefficients A^ and B^ 
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must satisfy the condition (A. 5) 

The magnetic potential in the cavity is 


-k'z 


- z + I C n e J 0 (k n p) 


where Jq is the Bessel function of order zero and k'a is 
the nth zero of jQ(k'a) = 0. From the gradient of (A. 35) 
we find the magnetic field 


= J, C n e k " Z k n J 0< k ,W 


’"1 k n C n e n J 0* k n o) 


Note that H p on the cavity boundary is identically zero. 

The magnetic field in the cavity can be written in terms 
of spherical coordinates using an identity from Cooke [1956] 
similar to the one employed in Appendix A 


IT £ C. f ( ■ 1 ) n ( k|r) n 
h; 1 - > ~ ) p (cose) + cose , 

t-1 n=l (n ‘ ,) - 


“ c Z (k'r) n dP o 

I T- Z(-D n 4-r ST- (cosej-slne 
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Equations (A. 38) and (A. 39) facilitate the application of the 
continuity conditions on the lunar dark side. 

Continuity of H r on r = a yields 


£ { -tA^a ' t/_1 - U+l P^( cos0 ) 

1 =] 


l ^ l M) 


(k;a) 


n v t 


1 = ] n = l 


(n-l)i P 


p„( cose ) 


COS0 = 


O<0< tt/2 , 


or 0 


TT 


< 0<TT 


From (A. 40) we find 


2 { l ' A i ,al _1 "(^ /+1 ) " 2 } ' f V,1 


( 2-t ' +1 ) 


(-D n (k'a) n 


V 1 r 

3 Z-J / 1 \ 

t-1 n=l (n - ,) 


l nl' 


^'=1 ,2,3 


whe re 


' nl ' (2l'+l ) 


if l ' =n 


dP 


dP./ 


V 3c5sff <°> ' P n BcosT «» 


(A. 40) 


(A. 41) 


n( n+T )-l '(-t'+l ) 


(A. 42) 
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Continuity of H 0 on r = a for O<0< tt/2 yields 

X (A t a t " , +B^a" t ' 2 )p[(cos9)=sin8+ f 1 ( k i a > np n ( cose > 

t=\ 1=] n=l 

, (A. 43) 


which upon expansion in terms of the set of odd associated 
Legendre functions can be written 


oo oo 

C' 


I ( A /’ 1+ V’ l ' 2 ) P U'"J 6 tM + I /I ('\ a)Up nt' ’ 

1=] 1=] n=l 


^'=1 ,3,5, 


(A. 44) 


Equations (A. 5), (A. 41) and (A. 44) determine A^, and 
for the aligned field case. In view of the detail presented 
in Appendix A, further algebraic reduction of the aligned 
field case will not be presented. 

The SSP sol ut ion is 


A ! ’ B 1 = 


aV 


, A^, B^=0 for 1*2, 


(A. 45) 


U I _ a^cose 

Hr OTT 


0 - 4 ) 

r 


(A. 46) 
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... -arsine ( , 

Hq " “TTt: V 1 


'0 


(a J -b J ) 


The SSV sol ution i s 



A,=l , Bj- b 3 / 2 , A^=B^=0 for 1*2 



(A. 47) 


(A. 48) 
(A. 49) 
(A. 50) 
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